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ABSTRACT 


This  report  documents  ocean  bottom  pressure  data  collected  from 


September  1983  to  May  1985  in  eleven  deployments  of  pressure  sensors 
under  the  Gulf  Stream  northeast  of  Cape  Hatteras  in  depths  of  3300  to 
3900  m,  as  part  of  the  Gulf  Stream  Dynamics  Experiment. 

In  past  experiments,  pressure  sensors  suitable  for  ocean  depths 
have  typically  exhibited  systematic  drifts  in  calibration  that  seriously 
contaminate  any  observed  periodicities  longer  than  a  few  days.  We  used 
Dlglquartz  sensors  (manufactured  by  Parosclentific,  Inc.),  because  these 
sensors  offered  potentially  much  lower  drift  than  other  commercially 
available  sensors.  In  these  sensors,  either  a  bellows  or  a  Bourdon  tube 
applies  stress  to  an  oscillating  quartz-crystal  beam,  causing  its 
oscillation  frequency  to  vary. 

Several  factors  Influence  the  amount  of  drift:  bellows  vs. 
Bourdon-tube  construction,  the  applied  pressure,  the  duration  of 
deployment,  and,  for  some  sensors,  high-pressure  preconditioning  in  the 
lab.  For  the  sensors  deployed  in  the  Gulf  Stream,  the  total  drift 
during  deployments  lasting  from  3  to  12  months  ranged  from  undetectable 
(So. 01  dbar)  to  0.20  to  0.50  dbar.  About  half  of  the  total  drift 
typically  occurred  within  the  first  6  days  of  deployment. 

We  estimate  the  residual  error  in  the  final  pressure  records,  after 
the  "dedrlftlng"  calculations,  to  be  typically  0.02  dbar  r.m.s.  (or 
0.06  dbar  r.m.s.)  if  the  first  6  days  of  the  record  are  excluded  (or 
included,  respectively).  This  low  drift-error  opens  many  possibilities 
for  studies  that  require  knowledge  of  the  low-frequency  dynamic  pressure 
signal  in  the  deep  ocean. 


Part  I  on  Methods  contains  a  short  review  of  bottom  pressure 


measurement  in  the  deep  ocean,  a  description  of  the  sensors  that  we 
used,  a  discussion  of  their  performance  and  drift  relative  to  type  of 
construction  and  prior  pressurization  history  ("preconditioning"),  and 
estimates  of  the  accuracy  of  the  dedrifted  pressure  records. 

In  Part  II  of  this  report,  the  full  data  processing  is  described, 
including  calibration  parameters,  corrections  for  the  influence  of 
temperature  variations  on  the  pressure  sensor,  and  parameterization  to 
remove  sensor  drift  errors  by  least-squares  regression  onto  an 
exponentially  decaying  time-dependence.  Time  series  are  plotted  which 
illustrate  several  steps  in  the  processing:  the  edited  half-hourly 
pressure  records,  the  detided  pressures  with  drift-model  curves 
superimposed,  and  the  low-pass  filtered,  "dedrifted"  pressure  records 
(i.e.,  after  subtracting  the  estimated  drift  curve). 
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Physical  oceanographers  have  long  desired  to  have  instrumentation 
capable  of  accurately  monitoring  the  small  variations  in  bottom 
pressures  in  the  deep  ocean  (Baker,  1981)  in  order  to  observe  the 
oceanic  pressure  field  in  analogy  with  the  use  of  barometers  by 
meteorologists.  The  measurement  problem  has  been  challenging  because 
the  dynamic  pressure  fluctuations  are  only  -0.10  to  0.01  dbar  in  an 
ambient  pressure  of  4000  dbar,  requiring  resolution  and  stability  of 
-1  ppm. 

In  the  last  fifteen  years,  deep-ocean  bottom  pressure  measurements 
have  been  conducted  in  several  studies  on  ocean  tides  (Filloux,  1971; 
Zetler  et  al.,  1975),  the  ocean  eddy  pressure  field  (Brown  et  al., 

1975),  weather-induced  bottom  pressure  fluctuations  (Beardsley  et  al., 
1977),  and  the  performance  of  the  pressure  instruments  themselves 
(Snodgrass  et  al.,  1975).  The  most  common  pressure  instruments  used  in 
the  above  studies  have  been  strain  gauges  (Wunsch  and  Wimbush,  1977), 
Vibrotron  sensors  (Wimbush,  1977),  metal  Bourdon  tubes  (Filloux,  1971; 
Mofjeld  and  Wimbush,  1977),  and  quartz-crystal  transducers  (Snodgrass  et 
al.,  1975). 

Although  the  requirement  of  adequate  sensitivity  and  short-term 

stability  is  of  great  importance,  the  chief  difficulty  in  dealing  with 

pressure  instruments  is  the  Instrumental  noise  of  low  frequency,  i.e., 

drift,  whose  magnitude  typically  increases  with  applied  pressure. 

Consequently,  it  has  been  most  difficult  to  measure  bottom  pressures  of 

the  deep  ocean.  The  drift  contaminates  the  observed  spectrum  most  at 

low  frequencies.  (A  simple  example  is  that  a  linear  drift  would  add  a 

-2 

red  noise  spectrum  of  slope  f  . )  Consequently,  earlier  studies  with 


deep  pressure  gauges  have  tended  to  focus  on  the  higher  frequency  (less 


contaminated)  aspects  of  the  data. 

In  Part  I,  we  present  the  experimental  setting  and  measurements 
(Section  2),  and  we  describe  the  pressure  and  temperature  sensors  (Section 
3).  We  next  discuss  the  observed  pressure  drifts  in  relation  to  sensor 
construction  and  to  prior  history  of  pressurization  (Section  4).  Next,  we 
describe  our  techniques  to  model  and  remove  the  drift  (Section  5),  and 
finally  we  attempt  to  estimate  the  uncertainties  in  the  drift  removal, 
i.e.,  the  accuracy  of  the  final  records  (Section  6).  All  the  time  series 
plots  of  the  data  are  shown  in  Part  II. 

2.  EXPERIMENTAL  SETTING  AND  MEASUREMENTS 

2.1  Experiment  Description  and  Overview 

As  part  of  the  Gulf  Stream  Dynamics  Experiment,  data  were  collected 
using  deep-ocean  bottom  pressure  and  temperature  sensors  in  the  Gulf 
Stream  northeast  of  Cape  Hatteras  from  September  1983  to  May  1985.  The 
measurements  were  made  under  the  combined  support  of  an  NSF  project 
entitled  "The  Dynamics  of  Gulf  Stream  Meanders"  and  an  ONR  project 
entitled  "Observations  on  the  Current  Structure  and  Energetics  of  Gulf 
Stream  Fluctuations  Downstream  of  Cape  Hatteras".  other  data  collected  as 
part  of  this  program  are  documented  in  separate  reports  (Tracey  et  al., 
1985;  Tracey  and  Watts,  1986b:  Frledlander  et  al.,  1986). 

The  deployment  sites  and  the  overall  study  area  are  shown  in  Figure 
2.1.  Inverted  echo  sounders  were  deployed  at  all  the  sites,  whereas 
pressure  sensors,  along  with  temperature  sensors,  were  deployed  only  at 
the  sites  indicated  by  the  solid  circles.  The  data  records  in  this  report 
are  labelled  by  (a)  the  letters  "PIES",  indicating  that  a  pressure  sensor, 
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a  temperature  sensor,  and  an  inverted  echo  sounder  were  deployed  at  the 
site,  (b)  the  last  two  digits  of  the  year  (84  or  85)  in  which  the 
instrument  was  recovered,  (c)  the  site  designation  indicated  in  Figure 
2.1,  and  (d)  the  letters  "CM"  between  the  line  letter  and  the  site  number 
if  a  current  meter  mooring  was  also  at  the  site.  For  example,  PIES84BCM2 
means  the  data  record  collected  in  1984  at  site  B2  where  a  current  meter 
mooring  had  also  been  deployed. 

The  time  series  plots  of  the  raw,  half-hourly  bottom  pressure  and 
temperature  measurements  are  given.  Processing  of  the  bottom  pressure 
records  Included  removal  of  the  tides  and,  in  particular,  the  long-term 
drift:  these  records  are  also  shown.  This  report  details  the  methods  used 
to  estimate  and  remove  the  drift,  and  to  estimate  the  residual 
uncertainty.  Subsequently,  both  the  bottom  pressure  and  temperature 
measurements  were  low-pass  filtered,  subsampled  at  six-hour  intervals,  and 
plotted.  First  order  statistics  (minimum,  maximum,  mean,  and  standard 
deviation)  are  listed  for  all  data  shown. 

_ Data  Recovery 

In  the  20-month  period  from  September  1983  to  May  1985,  seven 
different  pressure  instruments  (Table  2.1)  were  used  for  a  total  of  eleven 
deployments  at  six  sites  (Table  2.2).  These  instruments  were  launched, 
recovered,  and  redeployed  on  five  cruises  aboard  the  R/V  ENDEAVOR  (EN106, 
22-30  September  1983;  EN107,  1-3  November  1983;  EN118,  1-18  June  1984; 
EN124,  11-20  January  1985;  EN130,  7-21  May  1985)  and  one  cruise  aboard  the 
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R/V  OCEANUS  (0C144,  9-19  January  1984). 

In  the  eleven  deployments  of  the  deep-ocean  pressure  and  temperature 
sensors,  all  Instruments  were  recovered,  giving  an  Instrument  recovery 
rate  of  100%.  Due  to  a  malfunction  of  the  electronic  board  controlling 
one  pressure  sensor,  no  pressure  measurements  were  obtained  for  that 
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TABLE  2.1  Pressure  sensors  used 


Sensor 

Serial 

Number 

Model 

Working 

Range 

(dbar ) 

Type 

1 

8181 

75K-002 

0-3450 

bellows 

2 

8180 

75000 

0-3450 

bellows 

3 

17848 

76KB-032 

0-4100 

Bourdon-tube 

4 

17849 

76KB-032 

0-4100 

Bourdon-tube 

5 

17911 

76KB-032 

0-4100 

Bourdon-tube 

6 

18426 

76KB-032 

0-4100 

Bourdon- tube 

7 

19327 

46K-032 

0-4100 

Bourdon-tube 
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instrument  during  one  deployment  period  (record  PIES85CCM1 ) .  During  the 
same  deployment  period,  another  pressure  sensor  malfunctioned;  large  steps 
(both  positive  and  negative)  were  observed  in  the  measured  pressures 
(record  PIES85BCM3).  All  nine  remaining  Instruments  performed 
successfully,  giving  an  82%  data  return  for  the  pressure  measurements. 
Complete  records  were  obtained  for  all  eleven  deployments  of  the 
temperature  sensors;  thus  the  return  rate  was  100%  for  those  data. 


3.  DESCRIPTION  OF  SENSORS 


3.1  Pressure  sensor  and  mounting 

All  the  sensors  we  used  have  been  manufactured  by  Paroscientif ic, 

Inc.  A  detailed  description  of  the  digital  pressure  transducers  we  used 
is  given  by  Paros  (1976)  and  by  Wearn  and  Larson  (1982).  Briefly,  the  k,ey 
sensing  element  in  the  pressure  transducer  is  an  oscillating,  beam-shaped 
quartz  crystal,  which  is  piezoelectrically  induced  to  vibrate  in  its 
lowest  resonant  flexural  mode  by  an  oscillator  circuit.  The  oscillation 
frequency  of  the  crystal  varies  with  the  stress  load  applied  to  it  through 
a  lever  arm  attached  to  either  a  bellows  or  Bourdon  tube  (Figure  3.1). 

The  bellows  or  Bourdon  tube  is  pressurized  at  the  full  ambient  pressure  of 
the  ocean  via  a  long,  thin  capillary  tube  filled  with  mineral  oil.  The 
measurement  of  fluid  pressure  is  made  by  counting  the  output  frequency  of 
the  oscillator  circuit. 

The  pressure  sensors  were  powered  and  controlled  by  Sea  Data  Corp. 
model  XP35  electronics  cards  installed  in  inverted  echo  sounders  (lESs). 
The  lES  circuitry  and  mooring  configuration  are  described  by  Chaplin  and 
Watts  (1984).  These  instruments  were  tethered  less  than  1  m  above  the 
ocean  floor,  in  order  to  have  minimal  vertical  motion  of  the  instrument  in 


Figure  3.1  Schematic  of  the  bellows-type 
Dlglquartz  pressure  sensor  (Parosclentlf ic, 
Inc.).  In  the  Bourdon-tube  type  of  sensor, 
the  bellows  is  replaced  by  a  Bourdon  tube. 
The  sensor  Is  connected  to  the  environment 
through  a  small  capillary  tube  filled  with 
mineral  oil. 


response  to  bottom  currents;  we  estimate  less  than  0.5  cm  vertical 


excursion  in  the  peak,  observed  deep  current  speeds  (-30  cm/s)  in  this 
area. 

Another  potential  concern  for  bottom  pressure  measurements  is  the 
stability  of  the  mounting  on  the  sea  floor.  Most  evidence  indicates  no 
tendency  in  these  deployments  for  the  anchors  to  sink  into  the  mud  or 
for  the  bottom  to  slump  downslope:  all  but  one  of  the  observed  pressure 
drifts  discussed  in  the  following  section  are  toward  lower  (rather  than 
higher)  pressures.  Moreover,  the  acoustic  travel  times  (t)  measured 
simultaneously  on  these  lESs  show  no  indication  of  depth  change,  although 
only  large  or  sudden  depth  changes  would  be  likely  to  show  up  in  the 
r  records. 

3.2  Temperature  sensor 

Because  temperature  also  affects  the  oscillation  frequency  of  the 
quartz  crystal  in  the  pressure  sensor,  it  is  essential  to  Independently 
measure  the  temperature  of  the  crystal.  A  thermistor  (Yellow  Springs 
International  Corp. ,  model  44032)  was  placed  as  close  as  possible  to  each 
sensor. 

The  thermistor  was  powered  and  controlled  by  Sea  Data  Corp.  model 
DC37  electronics  cards  installed  in  the  lESs.  This  circuitry  produces  an 
output  frequency  that  varies  with  the  temperature-dependent  resistance  of 
the  thermistor. 

_ Time  base  and  recording 

The  frequency-counting  periods  for  the  above  measurements  were 
controlled  by  a  quartz-crystal  clock  in  the  lES;  the  circuitry  is 
described  by  Chaplin  and  Watts  (1984).  The  stability  of  the  crystal 
(JAN  crystals,  model  HC33/U  AT  cut  4.194304  MHz)  and  timing  circuitry 
used  in  these  deployments  was  reported  by  the  company  to  be  25  ppm/°C  and 


20  ppm  In  the  first  year.  However,  these  stability  specifications  are  not 
as  good  as  they  need  to  be;  the  JAN  crystals  were  mistakenly  used  as  a 
less  expensive  replacement  for  an  earlier,  better  crystal.  The  JAN 
crystals  had  all  been  baked  at  elevated  temperatures  and  aged  before 
deployment:  however,  clock-frequency  drift  could  still  be  responsible  for 
most  of  the  observed  drifts  in  the  measured  pressures.  (After  these 
deployments,  we  discovered  our  error  and  upgraded  the  time  bases.) 

All  the  data  were  recorded  digitally  on  the  lES  data  cassettes  using 
a  Sea  Data  model  610  recorder. 

4.  OBSERVED  DRIFTS  OF  PRESSURE  SENSORS 

lU _ General  description  of  performance 

In  many  records  a  long-term  drift  is  apparent  when  we  remove  the 
mean  and  the  tides.  Some  characteristics  of  this  drift  are  summarized 
in  Table  4.1,  and  plots  of  the  detlded  pressures  are  shown  in  Part  II 
(Figures  11,1  to  11.10). 

Several  features  are  common  to  these  records;  1)  During  the  first 
few  hours  (il2  h)  there  is  a  rapid  slewing  toward  higher  pressure  readings 
as  the  sensors  come  to  thermal  equilibrium  at  the  ocean  bottom.  2)  During 
the  next  3  to  12  months  of  deployment,  several  (five  of  the  ten)  pressure 
records  drift  by  several  tenths  of  a  decibar,  with  the  rate  of  drift 
decreasing  greatly  with  increasing  time.  3)  Superimposed  on  the  drift  is 
the  "ocean  pressure  signal",  with  obvious  5-  to  30-day  periodicities  and 
0.1-  to  0.2-dbar  pressure  changes;  these  are  the  dynamic  pressure  signals 
of  central  scientific  Interest  that  we  wish  to  study  uncontaminated  by 


drift.  4)  There  is  high-frequency  variability  of  roughly  ±0.02  dbar. 


Sensor  Preconditioni 


which  exceeds  measurement  noise  or  our  estimates  of  instrument  depth 


variation. 

Table  4.1  also  summarizes  the  drift  characteristics  found  for  all  the 
deployments;  for  simplicity,  only  the  total  change  during  the  first  6  days 
and  the  total  change  from  start  to  end  of  the  record  are  listed.  Several 
other  factors  are  also  listed  that  could  be  expected  to  Influence  the 
amount  of  drift.  These  include  the  ambient  pressure  (fairly  similar  at 
3340  to  3990  dbar  for  all  these  deployments),  the  sensor  construction 
(bellows  vs.  Bourdon  tube),  and  whether  the  sensor  was  "preconditioned"  to 
high  pressure  for  some  period  before  deployment.  The  following 
subsections  discuss  the  amount  of  drift  in  relation  to  these  factors, 
although  we  must  caution  that  we  have  very  few  replicate  cases  for 
comparison. 

4.2  Bellows  vs.  Bourdon-tube  constriictlon 

We  first  categorize  the  drift  performance  by  bellows  vs. 

Bourdon-tube  construction  (Table  4.1).  There  were  two  bellows  sensors 
with  three  deployment  records,  all  of  which  showed  significant  drift 
(-0.20  to  0.50  dbar).  Five  Bourdon-tube  sensors  used  in  eight  deployments 
produced  seven  records;  in  the  second  deployment  of  sensor  3  (record 
PIES85CCM1).  the  electronics  did  not  work.  In  four  of  the  seven  recovered 
records,  there  was  no  detectable  drift,  but  in  the  other  three  records, 
the  sensors  unfortunately  performed  unpredlctably :  One  sensor  (4;  record 
PIES84CCM2)  had  one  of  the  worst  drifts  of  all  (0.45  dbar)  in  its  first 
deployment  but  none  in  its  next  (record  PIES85C1).  (This  sensor  was  used 
in  two  different  lESs.  Clock-frequency  drift  could  be  responsible  for 
this  difference.)  Another  sensor  (7;  record  PIES85CCM3)  had  -0.20  dbar 
drift.  A  third  sensor  (6;  record  PIES85BCM3)  recorded  numerous 
discontinuous  Jumps  (-0.20  dbar)  of  varying  size  and  sign  (again,  possibly 


I 

i 
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attributable  to  the  time  base  crystal),  as  shown  in  Figure  11.3.  Only  one 
Bourdon-tube  sensor  (5)  did  not  drift  in  both  of  its  deployments. 

4.3  Effect  of  prior  pressurization  ("preconditioning") 

Secondly,  we  categorize  performance  according  to  whether  the 
sensors  were  "preconditioned"  at  high  pressure.  Both  bellows  sensors  were 
preconditioned.  One  bellows  sensor  (1),  which  was  pressurized  at 
3200  dbar  for  24  months  prior  to  its  deployment  (at  3625  dbar),  drifted 
only  about  half  as  much  as  the  other  bellows  sensor  (2).  That  sensor  was 
pressurized  at  3900  dbar  for  4  months  prior  to  deployment  (at  3730  dbar) 
and  drifted  with  opposite  sign.  (This  is  surprising,  because  Wearn  and 
Larson  (1982)  suggest  that  all  Internally  pressurized  bellows  sensors 
ought  to  drift  in  the  same  direction.)  The  opposite  direction  of  drift 
may  possibly  be  explained  by  the  fact  that  the  preconditioning  was  at  a 
slightly  higher  pressure  than  that  of  the  deployment. 

The  only  Bourdon-tube  sensor  (3)  that  was  preconditioned  showed  no 
drift,  but  other  Bourdon-tube  sensors  had  no  drift  without 
preconditioning. 


5.  METHOD  OF  DRIFT  ESTIMATION  AND  REMOVAL 

5.1  Introduction  to  drift  modeling 

Historically,  investigators  have  modeled  the  drift  of  various 
pressure  sensors  by  either  a  power-law  [-(t  -  t,)^,  0<^<1  ], 
logarithmic  ["log-law",  -ln(t  -  t^)  ]  or  exponential  {-exp[-a(t  -  t,)]} 
dependence  on  time,  t,  after  initial  pressurization  at  t  =  to.  Although 
the  exact  cause  of  the  drift  is  unknown,  it  is  believed  to  result  from 


mechanical  creep  of  materials  subjected  to  high  stress.  The  above  three 


dependencies  have  been  used  to  describe  creep  in  laboratory  and 
geophysical  studies,  as  has  been  reviewed  by  Wunsch  and  Wimbush  (1977)  and 
Wearn  and  Larson  (1982). 

It  had  been  hoped  and  asserted  (Wearn  and  Larson,  1980)  that  by 
careful  measurement  of  drift  characteristics  of  a  sensor  in  the 
laboratory,  the  drift  curve  for  it  in  an  open  ocean  deployment  might  be 
accurately  predicted.  Skepticism  has  probably  always  been  appropriate 
about  applying  such  predictions,  because  the  drift  processes  may  not  be 
reversible  or  reproducible.  Moreover,  for  the  deep-ocean  pressure  sensors 
it  has  not  yet  been  practical  to  monitor  the  drift  in  the  lab  at  the 
required  -1-ppm  accuracy  for  long  time  periods.  In  one  of  the  best 
attempts  (Wearn  and  Larson,  1980),  the  required  level  of  absolute  pressure 
and  temperature  stability  was  not  achieved  even  by  their  high-quality 
calibration  instruments. 

In  these  deployments,  the  approach  to  drift  removal  has  been  to  make 
least-squares  fits  of  the  above  mathematical  models  to  the  records 
themselves.  We  used  the  non-linear  least-squares  regression  routine, 

P3R,  in  the  BMDP-79  package  of  computer  programs  (Dixon  and  Brown,  1979). 
Since  the  records  are  long  compared  to  the  time  scales  of  the 
deep-ocean  pressure  variability,  the  "ocean  pressure  signal"  will  tend  to 
average  out  in  estimating  the  drift.  We  tried  fitting  our  data  to  all 
three  of  the  above  drift  dependencies  and  to  an  exponential-linear 
dependency:  the  results  are  compared  in  this  section.  The  accuracy  of  the 


best-fit  drift-curve  models  is  discussed  in  Section  6. 


.2  Log-law  vs. 


We  first  made  a  comparison  between  the  power-law  and  the  log-law 


types  of  drift,  described  by  the  following  formulas,  respectively; 

Vift  '  f. 


^drlft 

where  t  is  the  time  in  hours,  t,  is  the  time  when  the  drift  starts 
(13  hours  before  the  first  data  point  used,  as  explained  in  Section  8.3 
in  Part  II),  and  P, ,  P, ,  and  Pj  are  free  parameters  whose  values  are 
individually  determined  for  each  record  by  the  nonlinear  regression 
subroutine. 

The  residual  deviations  of  the  measurements  from  the  fitted  curves 
are  essentially  the  same  for  both  laws  and  are  shown  in  Table  5.1.  The 
two  types  of  curves,  when  plotted,  are  virtually  indistinguishable 
within  the  measurement  Interval.  The  log  law  is  preferred  because  of  its 
computational  simplicity  and  its  rapid  convergence  to  unique  parameters. 
(By  contrast,  the  regression  for  the  power  law  converges  slowly  and  the 
final  parameters  are  affected  by  the  starting  values  in  the  iteration 
process.)  Values  of  P,  and  P,  for  the  logarithmic  drift  curve  model  are 
listed  in  Table  5.2,  with  the  r.m.s.  deviations  repeated  for  each  record. 

The  logarithmic  curves  superimposed  on  our  drift  data  are  shown  in 
the  left-hand  column  of  Figure  5.1.  The  fitted  log-curves  lie  within  the 
scatter  of  the  data  and  approximate  the  existing  drift  remarlcably  well. 
However,  in  the  first  part  (-30  days)  of  the  records  PIES84CCM2  and 
PIES85BCM2,  the  predicted  drift  rates  are  probably  unreal  (overestimated 
for  PIES84CCM2  and  underestimated  for  PIES85BCM2). 


rABLE  5.1  The  r.m.s.  deviations  of  drift  curves  from  the  detided  pressures. 

Ihe  r..m.s.  AP .  =  [Z(AP  )^N  ,  where  AP.  is  the  difference  between  P,  . 

1  i  1  drift 

and  the  detided  pressure.  The  summation  extends  over  all  N  points  of  the 
tim.e  series.  The  r.m.s.  AP^  represents  the  r.m.s.  ocean  pressure  signal. 
Time  is  in  hours:  t.  is  the  time  when  the  drift  starts. 
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r.m.s.  AP^ 

(dbar ) 

Exponential- 

Power 

Logarithmic  Exponential  Linear 

Record 

Sensor 

Law 

Law  Law  Law 

PIES84B2 

PIES85BCM2 

PIES85CCM2 

PIES84CCM2 

PIES35CCM3 


0.0372 

0.0491 

0.0361 

0.0543 

0.0412 


0.0428 

0.0474 

0.0363 

0.0542 

0.0416 


0.0372 

0.0368 

0.0421 

0.0448 

0.0370 


0.0370 

0.0368 

7.0349 

0.0449 

0.0370 


TABLE  5.2  Drift  parameters  for  log-law  least-squares  fit: 


'■drift  = 

or,  in  the  case  of  no  drift,  Pj  is  zero. 

drif  t 

Parameters  as  in  Table  5.1. 


Record 

Sensor 

Pi 

(dbar ) 

Pi 

( dbar ) 

r.m.s.  AP^ 

(dbar ) 

PIES84B2 

1 

-0.037278 

0.231444 

0.043 

PIES85BCM2 

1 

-0.048840 

0.394873 

0.047 

PIES85CCM2 

2 

0.088609 

-0.709164 

0.036 

PIES84C1 

3 

— 

— 

0.041 

PIES84CCM2 

4 

-0.112501 

0.852820 

0.054 

PIES85C1 

4 

— 

— 

0.030 

PIES84CCM3 

5 

— 

— 

0.053 

PIES85C0 

5 

— 

— 

0.042 

PIES85CCM3 

7 

-0.035110 

0.281740 

0.042 

5.3  Exponential  and  exponential-linear  vs.  log-law  least-sauares  fit 

We  next  made  least-squares  fits  of  exponential  and  exponentlal-linear 
curves  to  our  data.  The  mathematical  formulas  we  used  were: 

^drift"  ~ 

for  the  exponential  curve  and 

Pdrift"  ■  exp(P,  t)]  +  P,  +  P,  t 

for  the  exponential-linear  curve.  Here  P, ,  P, ,  P, ,  and  P4  are  again  free 
parameters  determined  by  the  non-linear  regression  subroutine  and  t  is 
time  in  hours,  relative  to  13  hours  before  the  first  data  point  used  (as 
explained  in  Section  8.3  of  Part  II).  (N.B.:  Both  models  could  use  any 
other  chosen  start  time,  to,  since  any  time  shift,  t  -  t,,  could  be 
adjusted  by  corresponding  changes  in  P,  and  Pj.) 

The  above  models  converge  rapidly  and  the  parameters  are  independent 
of  initial  values  in  the  iteration  process.  Values  of  P, ,  P, ,  P,  and  P, 
for  the  exponential  and  exponential-linear  drift  curve  models  are  listed 
in  Table  5.3,  with  r.m.s.  deviations  repeated  for  each  record.  The 
predicted  exponential  and  exponential-linear  drift  curves  have  in  general 
smaller  r.m.s.  deviations  from  the  original  data  than  the  logarithmic 
curves  (Table  5.1).  This  was  expected  since  they  have  more  free 
parameters  (four).  However,  curves  for  all  the  laws  have  similar  shapes, 
and  the  r.m.s.  residual  is  mainly  from  the  oscillatory  "ocean  pressure 
signal".  Plots  of  exponential  and  exponential-linear  curves  superimposed 
on  the  drift  data  are  shown  in  the  center  and  right-hand  columns  of 
Figure  5.1,  respectively. 

The  regression  subroutine  detects  no  linear  trends  in  records 
PIES85BCM2,  PIES84B2  and  PIES85CCM3.  It  gives  P,  =  0  and  converges  to  the 
same  values  for  P, ,  P, ,  and  Pj  with  the  purely  exponential  type  of  drift. 
For  these  three  records  the  exponential  curves  seem  very  satisfactory  in 
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traclcing  the  drift  and  were  preferred  over  the  logarithmic  ones.  In 
records  PIES84CCM2  and  PIES85CCM2  we  observe  a  probable  tendency  for  the 
exponential  law  to  underest i.mate  the  drift  rates  at  the  beginning  of  the 
records  and  for  the  exponential-linear  law  to  overestimate  the  drift  rates 
at  the  end  of  the  records.  In  the  case  of  PIES84CCM2  we  preferred  the 
exponential  curve:  with  its  fewer  parameters,  the  r.m.s.  deviation  was  the 
best:  we  chose  to  accept  the  risk,  of  retaining  a  small  additional  drift 
error  in  the  first  -10  days  of  the  record.  In  the  case  of  PIES85CCM2  we 
preferred  the  exponential-linear  curve;  the  r.m.s.  deviation  was 
significantly  smaller  with  this  law  and  we  chose  to  accept  the  risk,  of 
introducing  a  linear  drift  trend  of  unknown  origin  (perhaps  due  to  the  JAN 
crystals ) . 

5.4  The  final  ("residual")  pressure  records 

The  "residual"  pressure  records  (detided  and  dedrifted),  which  are 
our  best  estimate  of  the  true  deep-ocean  pressure  signal,  were  calculated 
by  subtracting  the  fitted  curve,  represented  by  the  fit  described 

above  (Table  5.3),  from  its  respective  detided  set  of  measurements  at  each 
half-hourly  data  point. 

A  low-pass  filtered  version  of  the  pressure  records  is  given  in 
Figures  12.4  to  12.6  of  Part  II.  All  coincident  records  are  plotted  on 
the  same  time  axes  and  with  consistent  pressure  scales.  (January-May  1984 
is  plotted  twice  to  show  the  coincidence  with  earlier  and  later  deployment 
periods.)  Visually,  the  coherence  between  pressure  records  at  neighboring 
sites  is  quite  high.  We  have  Identified  some  of  the  major  pressure-events 
with  the  passage  of  Gulf  Stream  rings  offshore  and  their  coalescence  into 


the  Gulf  Stream  (Bane  and  watts,  1985). 


6.  ESTIMATED  ACCURACY  OF  FINAL  PRESSURE  RECORDS 


6.1  General  remarks 

The  final  record  described  in  the  previous  section  is  a  pressure 
record  from  which  a  drift  curve  has  been  subtracted;  the  best  estimated 
drift  curve  is  the  exponential  or  exponential-linear  least-squares  fit 
whose  parameters  are  listed  in  Table  5.3.  We  now  attempt  to  estimate  the 
uncertainty  in  the  drift-removal  procedure  or,  equivalently,  to  estimate 
the  long-term  accuracy  of  the  dedrifted  records.  (N.B.:  This  is  a 
separate  question  from  that  of  the  absolute  accuracy  of  the  mean  pressure, 
which  is  about  ±1.5  dbar.  The  absolute  accuracy  has  not  been  important 
thus  far,  because  the  absolute  depth  of  the  instrument  site  relative  to 
the  geold  is  not  Rnown  either.) 

It  is  difficult  to  estimate  the  accuracy  of  the  drift  curves  for 
these  records.  Each  has  a  mixture  of  two  unknown  signals  that  differ  for 
each  site:  the  ocean  pressure  signal  plus  the  drift.  Some  records  are 
from  different  time  periods,  and  others  coincide  in  time  but  are  from 
different  sites  (-15  to  50  km  apart)  with  coherent,  but  different  ocean 
pressure  signals.  We  have  no  independent  way  to  determine  either  the 
ocean  pressure  signal  or  the  drift  signal  exactly,  but  the  methods  we  have 
used  to  separate  them  essentially  assume  that  they  have  different  spectral 
characteristics . 

Since  the  rate  of  drift  decreases  greatly  with  increasing  time, 
another  distinction  is  that  the  longer  duration  records  (especially  their 
later  months)  have  a  smaller  proportion  of  their  total  variance  associated 
with  drift  errors.  Hence,  they  become  more  dominated  by  the  ocean 


pressure  signal. 


We  take  advantage  of  these  differences  to  estimate  the  errors  in  two 
ways.  In  the  first  method,  we  use  records  from  two  sites  at  which  an 
instrument  that  had  been  deployed  for  several  months  was  recovered  and 
replaced  almost  immediately  with  another  instrument  at  the  same  site. 
Figure  12.5  in  Part  II  shows  these  cases  (records  PIES84CCM2  and 
PIES85CCM2  at  site  C2;  records  PIES84CCM3  and  PIES85CCM3  at  site  C3),  with 
the  recovery  and  redeployment  time  in  June  1984.  The  drift-removal  error 
at  the  end  of  one  record  is  probably  small  compared  to  that  at  the 
beginning  of  the  subsequent  record,  since  it  seems  reasonable  to  expect 
that  the  magnitude  of  error  is  roughly  proportional  to  the  amount  of  drift 
being  removed.  Hence,  any  pressure  "Jump"  occurring  between  the  records 
gives  an  estimate  of  the  size  of  drift-removal  error  that  characterizes 
the  beginning  of  a  record.  The  "jump"  magnitudes  are  about  0.02  dbar  at 
site  C2  and  0.06  dbar  at  site  C3  (also  listed  in  Table  6.1A). 

The  second  method  to  estimate  the  error  was  to  regress  the  law 
selected  for  dedrifting  on  segments  of  the  records.  These  segments  have 
the  same  (true)  drift  curve,  but  our  estimate  of  it  will  differ  for 
different  segments  due  to  the  contribution  of  the  ocean  pressure  signal. 

An  inferrence  of  how  much  the  ocean  pressure  signal  "contaminates"  or 
causes  the  regressed  drift  curve  to  differ  from  the  true  drift  can  be 
taken  from  the  difference  between  the  regressed  drift  curves  obtained  from 
different  segments  of  the  record. 

We  used  the  second  method  only  on  the  four  longest  records  (7  to  12 
months),  because  it  was  felt  that  dividing  the  records  into  segments 
shorter  than  3  months  was  inadvisable,  given  the  obviously  energetic  ocean 
pressure  signals  at  shorter  periodicities.  Each  of  the  four  records  was 
divided  into  two  segments  at  a  point  which  had  a  pressure  near  the  overall 


TABLE  6.1  Pressure  offsets  between  drift  curves  estimated  on  different 

sequential  records  (6.1A)  and  on  segments  of  the  same  record  (6. IB); 

estimated  accuracy  of  final,  "dedrifted"  pressure  records.  The 

detection  threshold  (due  to  "noise"  of  the  ocean  pressure  signal) 

of  the  exponential  and  exponential-linear  least-squares  fit  is 

<  0.01  dbar  in  1  year.  AP  =  AP.  -  AP  . 

b  a 

t  :  Length  of  first  segment  or  record 
t  :  Length  of  second  segment  or  record 

^P:  Pressure  offset  between  drift  curves  of  two  segments  at  Join-point 

or  of  two  full  records 

AP  :  Pressure  offset  between  drift  curve  of  first  segment  and  drift 
^  curve  of  the  full  record  at  the  end  of  first  segment 
AP  :  Pressure  offset  between  drift  curve  of  second  segment  and  drift 
curve  of  the  full  record  at  the  beginning  of  second  segment 


6  .  lA 

Pressure  offsets  between 

different 

,  sequential 

records 

t 

a 

^b 

AP 

Records 

(months ) 

(months ) 

(dbar ) 

PIES84CCM2 

and 

8 

11 

0.02 

PIES85CCM2 

PIES84CCM3 

and 

6 

11 

0.06 

PIES85CCM3 

6.  IB 

Pressure  offsets  between 

segments 

Of  the  same 

record 

t 

a 

’^b 

AP 

AP 

a 

Record 

(months ) 

(months) 

( dbar ) 

(dbar ) 

(dbar ) 

PIES85BCM2 

4.4 

8 

-0.015 

0.013 

-0.002 

PIES85BCM2 

7.4 

5 

0.000 

0.001 

0.001 

PIES84CCM2 

3 

4.3 

-0.016 

-0.052 

-0.036 

PIES84CCM2 

4.6 

2.6 

0.042 

-0.020 

0.022 

PIES85CCM2 

2.6 

8.6 

-0.007 

-0.004 

-0.011 

PIES85CCM2 

5.6 

5.6 

0.037 

-0.028 

0.009 

PIES85CCM2 

8.3 

2.9 

0.053 

-0.004 

0.049 

PIES85CCM3 

2.6 

8.6 

0.003 

0.014 

0.017 

PIES85CCM3 

5.6 

5.6 

-0.006 

0.006 

0.000 

r . m . s .  AP 

for 

segments : 

0.027 

0.021 

0.022 

r . m . s .  AP 

for 

full 

12-month 

record* : 

0.019 

0.015 

0.016 

« / 

*  Obtained  by  dividing  by  (2)  *.  because  the  degrees  of  freedom  are 

doubled. 


drift  curve,  to  avoid  introducing  new  "end  effects".  A  curve  was  fitted 
to  model  the  drift  for  each  segment.  Then  we  computed  the  offsets  between 
these  fitted  curves  at  the  point  where  one  Joins  or  overlaps  another.  At 
the  same  point,  we  also  computed  the  differences  between  the  curves  fitted 
on  the  segments  and  the  curve  originally  fitted  on  the  whole  record.  The 
results  of  these  calculations  are  summarized  in  Table  6. IB.  The 
magnitudes  of  the  offsets  range  from  0.000  to  0.053  dbar,  with  a  r.m.s. 
offset  of  0.023  dbar. 

6.3  Discussion  of  drift  uncertainty  estimates 

We  now  argue  that  the  offsets  between  the  above  exponential-law 
regressions  on  segments  of  each  given  record  are  probably  overestimates  of 
the  uncertainty  of  the  drift  regression  on  the  full  record.  Two  reasons 
support  this  contention:  (1)  The  offsets  occur  at  the  beginning  of  each 
segment,  where  the  regressed  curve  has  its  maximum  amplitude  and  probably 
maximum  error.  However,  in  the  regression  curve  for  the  full  record,  the 
corresponding  time  is  not  at  the  beginning,  but  farther  out  in  the  tail  of 
the  true  drift  curve,  where  the  drift  rate  is  slower  and  the  error  should 
be  smaller.  (2)  The  full-record  regression  curve  is  2  or  4  times  as  long 
as  the  Independent  segments,  with  consequently  greater  degrees  of  freedom 
and  more  confidence  in  "averaging  out"  the  ocean  pressure  signal;  invoking 
a  "central-limit  theorem"  that  the  ocean  signal  ought  to  be  normally 
distributed  given  sufficient  degrees  of  freedom,  would  suggest  that 
doubling  (or  quadrupling)  a  record  length  would  Improve  the  removal  of 
ocean  pressure  signal  from  the  drift  by  a  factor  of  2  ’  (or  2, 

respectively).  These  arguments  are  applied  in  the  footnote  to  Table  6. IB. 


7.  SUMMARY  OF  PART  I 


Ten  records  from  seven  Dlglquartz  deep-ocean  bottom  pressure  sensors 
( Parosclentif ic ,  Inc.)  have  been  obtained  in  eleven  deployments  of  3  to  12 
months  between  September  1983  and  May  1985.  The  deployments  were  under 
the  Gulf  Stream  northeast  of  Cape  Hatteras  in  depths  of  3300  to  3900  m,  as 
part  of  the  Gulf  Stream  Dynamics  Experiment.  Two  of  the  pressure  sensors 
had  bellows;  five  had  Bourdon  tubes.  The  records  were  examined  with 
particular  attention  to  estimating  and  removing  any  observed  drift  in 
calibration.  (For  this  set  of  deployments,  much  of  the  drift  may 
unfortunately  be  due  to  drift  of  the  control  crystals  for  the  time  base. 
Using  better  clock,  crystals  in  future  deployments  should  result  in  less 
drift.  ) 

All  three  records  from  the  two  bellows-type  sensors  exhibited 
significant  drift.  (Surprisingly,  one  sensor  drifted  in  the  opposite 
sense  of  the  other  sensor.) 

Seven  records  were  obtained  with  the  five  Bourdon-tube  sensors.  Two 
of  these  sensors  were  deployed  twice;  of  these  four  records,  three  records 
had  no  detectable  drift  and  one  record  (a  first  deployment)  had  a  large 
drift.  Three  sensors  were  deployed  once;  of  these  three  records,  one  had 
no  detectable  drift,  one  had  significant  drift,  and  one  was  Jumpy  and 
unreliable. 

The  drift  was  removed  from  five  pressure  records  as  follows: 
Least-squares  fits  of  the  power-law,  logarithmic,  exponential,  and 
exponential-linear  dependencies  on  time  were  made  to  each  detlded  pressure 
record.  The  best-fit  drift  curve  was  subtracted  from  the  record  at  each 
data  point,  yielding  the  "residual"  pressure  record.  For  four  records. 


the  exponential  drift  curve  had  the  best  fit;  for  one  record,  the 


exponential-linear  drift  curve  had  the  best  fit. 

We  estimate  that  the  errors  in  our  drift-removal  procedure  are  as 
much  as  0.06  dbar  (range  ±0.02  to  ±0.06  dbar;  Table  6.1A)  if  the  entire 
record  is  considered.  However,  the  pressure  error  is  estimated  to  be 
about  0.02  dbar  r.m.s.  (range  ±0.01  to  ±0.05  dbar;  Table  6. IB)  if  the 
first  week,  of  a  record  is  not  included. 

The  residual  uncertainties  in  the  final  "dedrlfted"  pressure 
records  are  small  compared  to  the  r.m.s.  ocean  pressure  signal  observed, 
and  will  tend  to  contaminate  only  the  relatively  long  periodicities 
(^60  to  90  d)  in  these  records.  This  opens  many  possibilities  for  studies 
requiring  knowledge  of  the  deep-ocean  dynamic  pressures.  Moreover,  bottom 
pressure  gauges  In  combination  with  Inverted  echo  sounders  can  provide 
highly  accurate  definition  of  variations  in  the  sea-surface  height. 
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8.  DATA  PROCESSING 


Flow  Chart  and  Description 


The  data  processing  for  all  the  lES,  pressure,  and  temperature 
records  was  done  on  a  PRIME  750  computer,  using  FORTRAN  77  with  some 
extensions  for  bit  manipulation.  The  basic  procedural  steps  on  the 
pressure  data  are  illustrated  by  the  flow  chart  in  Figure  8.1  and  are 
accomplished  by  a  series  of  routines  developed  specifically  for  the  lES 
(Tracey  and  Watts,  1987).  The  raw  data  are  recorded  within  the  lES  on 
Sea  Data  model  610  recorders.  The  cassette  tape  contains  the  counts 
associated  with  sequential  record  number,  acoustic  travel  time, 
pressure,  and  temperature  measurements  as  a  series  of  integer  words  of 
varying  lengths.  The  subsequent  processing  steps  are  outlined  below. 

CARP:  Transfers  the  data  from  cassettes  to  9-track,  magnetic  tape  for 
subsequent  processing. 

BUNS:  Converts  the  series  of  integer  words  of  varying  lengths  into 
standard  length  32-blt  integer  words. 

MEMOD:  Establishes  the  time  base  from  the  sequence  number.  Converts 
all  travel  time,  pressure  and  temperature  counts  into 
scientific  units  of  seconds,  decibars,  and  degrees  Celsius, 
respectively. 

FILL:  Checks  for  proper  increment  of  the  time  base.  Missing  data 
points  are  filled  by  inserting  interpolated  values. 

DESPIKE:  Identifies  and  replaces  travel  time,  pressure  and 

temperature  spikes  with  Interpolated  values. 

RESPO:  Determines  the  mean  and  the  tidal  constituents  by  using  tidal 
response  analysis  (Munk  and  Cartwright.  1977).  The  mean  and 


the  calculated  tides  are  then  removed  from  the  pressure 
values . 

DEDRIFT:  Identifies  and  removes  the  long-term  drift  of  the  detided 

pressures.  The  drift  is  tracked  through  the  non-linear 
regression  subroutine  P3R  of  BMDP-79,  a  package  of  computer 
programs  developed  at  UCLA  (Dixon  and  Brown,  1979). 

LOW-PASS  FILTERING:  Convolves  the  residual  (detided  and  dedrifted) 
pressures  with  a  40-hour  low-pass  Lanczos  filter.  The  FESTSA 
time  series  analysis  package  (Brooks,  1976),  modified  for  the 
PRIME  750,  was  used  to  remove  the  higher  frequency  (tidal  and 
Inertial)  motions  from  those  with  periods  of  several  days  or 
longer.  The  symmetric  filter,  with  a  Lanczos  taper,  was 
designed  with  the  half  power  (-6  db)  at  0.025  cph  and  -60  db 
attenuation  at  tidal  frequencies  (0.0833  cph).  The  smoothed 
output  time  series  ("40HRLP")  for  pressure  and  temperature  were 
subsampled  at  six-hour  intervals. 

8.2  Conversion  to  Scientific  Units 

Both  the  temperatures  (T)  and  the  pressures  (P)  are  determined  by 
counting  frequencies  that  vary  in  accord  with  individually  calibrated 
relationships  for  T  and  P.  These  are  described  below  and  the  individual 
parameters  for  each  Instrument  are  given  in  the  tables  in  Section  9. 
8.2.1  Temperature 


The  T  measurement  is  independent  of  P;  that  is,  the  thermistor  does 
not  vary  with  pressure.  Thus,  only  the  temperature-dependent 
resistance,  R(T),  of  the  thermistor  controls  a  variable-frequency 
oscillator  on  the  DC37  electronics  card  (Sea  Data  Corp.).  This 
circuitry  is  designed  such  that  it  linearizes  the  output  frequency,  f^, 
as  a  function  of  T  within  the  range  0°C  to  25”c.  (This  linearization  is 


5b 


set  for  a 


R(T)  curve  for  this  model  YSI  thermistor  and  is 


supposed  to  be  accurate  to  0.1  C. )  The  frequency,  f^,  is  counted  for  a 

measurement  interval,  Z,  and  divided  by  a  scale  factor,  S  (=  1  in  this 

case),  to  fit  into  a  16-bit  integer  word  prior  to  recording  on  the 

cassette  tape.  Therefore,  the  recorded  count  is  =  Zf  /S.  From  N 

T  T  T 

one  can  find  T  (°C),  which  is  the  true  average  temperature  over  the 
sampling  Interval,  Z.  For  these  deployments,  the  sampling  interval  was 
30  min  (1800  s). 


The  conversion  from  counts  to  temperature  is 

'“t  ■  -  "a>  ■  ^  -  ■'a>  ■ 

where  the  constants  and  (T  .r  ),  listed  In  Table  8.1  for  each 

a  b  a  b 

sensor,  are  determined  from  laboratory  calibrations.  The  accuracy  of 

the  temperature  measurements  is  about  0.1°C  and  the  resolution  is 
o 

0,0002  C,  The  half-hourly  temperature  records  are  shown  in  Section  10 


(Figures  10.1-10).  The  temperatures  were  low-pass  filtered  and  are 
shown  in  Section  12  (Figures  12.1-3). 


The  pressure  sensor,  suited  for  deep-ocean  applications,  has  an 
output  frequency,  fp,  that  not  only  depends  strongly  on  P,  but  also 
depends  significantly  on  T.  The  full-range  change  in  frequency  is 


typically  from  36  to  40  kHz,  depending  on  the  applied  pressure  and  type 
of  construction.  It  is  common  practice  to  deal  with  the  period, 


T  =  1/fp,  rather  than  the  frequency;  the  period  ranges  from  25  to 
28  Msec.  Accompanying  each  sensor,  we  purchased  individual  calibrations 
from  Paroscientif Ic ,  Inc.  along  with  polynomial  expressions  to  determine 
P(r,T). 

An  XP35  electronics  card  (Sea  Data  Corp. )  counts  the  pressure 
sensor  oscillations  for  the  measurement  Interval,  Z,  and  divides  the 


insignificant. 
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counts  by  the  scale  factor,  S,  to  fit  into  a  24-bit  integer  word  for 
recording  on  the  cassette  tape.  The  sampling  interval  was  30  min 
(1800  s)  and  the  scale  factor  was  4  for  these  deployments.  Therefore, 
the  recorded  count  is  Np  =  Zfp/S  =  450/r. 

The  polynomial  relationships  P(t,T)  from  Paroscientif ic ,  Inc.  are 
listed  in  Table  8.2  for  each  sensor,  t,  is  the  vibrational  period  at 
zero  applied  pressure,  with  polynomial  temperature  dependence.  The 
older  sensors  come  with  "A,B,t,"  coefficients,  each  of  which  has  its  own 
polynomial  equation  in  temperature,  T  (  F).  The  newer  sensors 
preferably  have  "C,D,ra"  coefficients,  again  each  with  a  polynomial  T 
dependence  (again  in  °F,  except  sensor  7,  which  is  in  °C). 

The  half-hourly  measured  pressure  records  are  shown  in  Figures 
10.11  to  10.20.  The  absolute  accuracy  of  the  long-term  average  pressure 
measurements  is  within  1.5  dbar,  but  the  resolution  is  much  better.  The 
pressure  change  corresponding  to  a  change  in  Np  of  one  count,  under 
typical  conditions  of  applied  pressures  (-3500  dbar)  and  deep 
temperatures  (-2°C),  ranges  from  0.0010  to  0.0025  dbar  for  these 
deployments  with  a  measurement  interval  of  Z  =  1800  s.  The  stability  of 
these  pressures  is  discussed  in  Section  8.3. 


The  date  and  time  were  assigned  to  each  sampling  period.  The 
tables  in  Section  9  report  the  hour,  minutes,  and  seconds  associated 
with  the  first  and  last  sampling  period  as  a  six-digit  number.  All 


times  are  given  as  Greenwich  Mean  Time  (GMT).  For  processing 
convenience,  the  times  were  converted  into  yearhours.  Table  8.3  lists 
the  yearhour  that  corresponds  to  0000  GMT  of  each  day  for  non-leap 
years.  (For  leap  years,  the  yearhours  can  be  determined  by  adding  24  to 
each  day  after  February  28.)  There  are  a  total  of  8760  hours  in  a 
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TABLE  8.2  Pressure  sensor  calibration  equations,  P(r.T),  to  convert 
sensor  oscillation  period,  r  (s),  into  pressure,  P  (psi),  with 
functional  dependence  on  temperature,  T  (either  ®F  or  ‘"C),  and  on 
the  output  period  at  zero  pressure,  r,  (s). 


Set  1 :  (for  sensors  1  and  2) 


P  =  A(l-  To/r)  -  B(l-  r./r] 


B  =  B, 


+  A,  T  +  A,  T* 
+  B,  T  +  B,  T’ 
+  r,  T  +  Tj  T* 


T  =  450/Np 

T  :  (°F)  (Fahrenheit) 


A„:  5.18004’10 
Ai :-9. 70308*10' 
Aj  :  1.71739*10' 

A,:  4.87536*10'^ 
Aj  :-9.17509*10' 
A,  :-5.94345*10' 


Sensor;  1.  Serial 
B,:  3.17505*10_ 
Bj:-7.80773*10_p 
Bj  :  1.04970*10  ^ 
Sensor;  2.  Serial 
Bo*.  2.81395*10^ 
B,;  8.01378 
B.  ;-9. 13550*10  ^ 


,;25. 97996*10 

-1.99543*10] 

1. 70393*10* 

,  -.25.95124*10' 
-6.38494*10' 
2.28777*10* 


Set  2 ;  (for  sensors  3,  4,  5,  6) 


P  =  C[(1-(to/t)*  )  -  D(l-(T,/r)M*] 


C  =  Co  +  Cl  T  +  Cj  T* 
D  =  Do  +  D,  T  +  Dj  T* 
ro=  Too  +  T,  T  +  r,  T* 


T  =  450/Np 

T  ;  (°F)  (Fahrenheit) 


Sensor; 


Co:-2. 95074*10 
C,  •.-1.40724*  10' 
C,;  6.61457*10* 


Co  :-2, 48035*10 
Cj;  5.17487*10* 
Cj  ;-4. 92268*10* 

Co  :-2.27081*10'’ 
C,  .--5.35951*10* 
C,  ;  2.09832*10* 


o;-2. 64209*10 
, ;-2. 91091*10] 
;  1.09997*10* 


D,:  5.00619*10 
Di:  4.19904*10* 
D, ;-l. 20873*10* 
snsor;  4.  Serlaj 
Do*.  4.55063*10 


Di: 

D,  ; 
Sensor 


1.06118*10 
-2.66707*10* 
:  5.  Serlal_ 
5.62671*10 
-3.31430*10* 


D, ;-3.31430*10_ 
D,;  1.16763*10* 
Sensor;  6  Serial 
D,;  6.22649*10_ 
D, :-2. 16853*10* 
D. ;  7.50013*10* 


;28. 05943*10 
-1.72125*10* 
6.73124*10* 

,  ;28.92787*10] 
-1.45426*10] 
7.89478*10* 


:29.75394*10_ 

-1.95993*10* 

7.66406*10* 


;29.16313*10_ 

-1.94124*10] 

5.35286*10* 


(for  sensor  7) 


P  =  C[(l-(T./r)*)  -  D(l-(r,/T)M*] 


C  =  Co  C,  T  +  C,  T*  T  =  450/N 

D  =  Do 

To  =  Too-*-  T,  T  +  T,  T*  -t  Tj  T’ ,  T  :  (°C)  (Centigrade) 


2.505732*10] 

6.270068*10* 

5.399068*10* 


Do:  3.720884*10 


Too :29.613417*10_ 
Tj;  -5.454051*10* 
T,;  3.093264*10] 
T,;  8.444666*10* 


yearhour  corresponding  to  0000  GHT  is  listed  for  each  day. 
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standard  year  and  8784  hours  in  a  leap  year.  The  yearhours  given  in 
this  report  are  referenced  to  0000  GMT  on  either  January  1,  1984  or 
January  1,  1985,  depending  on  the  year  in  which  the  instrument  was 
recovered:  the  two-digit  number  of  the  instrument  name  indicates  which 
date  is  the  reference.  Positive  yearhours  correspond  to  sampling 
periods  which  occur  during  the  same  calendar  year  as  the  reference  date; 
negative  yearhours  correspond  to  those  which  occur  in  the  calendar  year 
prior  to  the  reference  date. 

8.3  Pressure  Sensor  Drift  Estimation  and  Removal 

Prior  to  determining  the  long-term  drift  of  our  measurements,  the 
mean  and  the  tides  were  removed  from  the  original  pressure  measurements. 
The  time  origin  (when  the  drift  started;  t,,)  was  assumed  to  be  the  time 
when  the  instrument  was  about  halfway  to  the  sea  floor,  or  1  hour  before 
the  first  sample  on  the  bottom.  We  also  removed  the  first  12  hours  of 
data  (24  points)  after  the  Instrument  landed  on  the  sea  floor  because 
during  this  time  the  sensors  inside  the  glass  instrument  housing  were 
still  coming  to  thermal  equilibrium,  causing  these  initial  measurements 
to  change  rapidly.  Hence,  the  time  origin  was  13  hours  before  the  first 
data  point  used.  For  computational  ease  and  speed  in  estimating  the 
drift  curve,  we  temporarily  subsample'’  the  time  series  every  two  hours: 
however,  the  drift  curve  was  later  removed  from  the  half-hourly  time 
series . 

The  method  of  estimating  and  removing  the  long-term  drift  from  the 
pressure  records  was  described  in  Section  5,  and  the  fitted  drift 
parameters  were  listed  in  Table  5.3.  The  regressed  drift  curves  are 


superimposed  on  the  detided  pressures  In  Figures  11.1  to  11.10. 


The  residual  pressure  records  (detided  and  dedrifted),  representing 

our  best  estimate  of  the  true  deep-ocean  pressure  signal,  were 

calculated  as  the  result  of  subtracting  the  fitted  P,  curve  from  the 

drift 

detided  sets  of  measurements  at  each  half-hourly  data  point.  These 
dedrifted  pressures  were  then  low-pass  filtered,  and  are  shown  in 
Figures  12.4  to  12.6. 


9.  PARAMETERS  AND  STATISTICS  OF  INDIVIDUAL  SITES  AND  RECORDS 


The  location,  dates,  and  basic  statistics  of  the  eleven  instrument 

records  are  given  in  Tables  9.1  to  9.11.  Each  table  documents  a  single 

instrument  deployment.  The  useful  data  records  are  plotted  in  Sections 

10,  11,  and  12.  General  site  information,  such  as  position,  bottom 

depth,  and  launch  and  recovery  times,  are  given  first,  followed  by 

details  about  the  bottom  pressure  and  temperature  data  sets.  The  times 

associated  with  the  first  and  last  data  points  are  supplied.  The  time 

origin  (to)  is  13  hours  before  the  first  data  point  used.  Yearhours  are 

referenced  to  0000  GMT  on  January  1  of  the  year  indicated  by  the 

two-digit  number  in  the  instrument  name.  Positive  yearhours  occur 

during  the  same  calendar  year  as  the  two-digit  number;  negative 

yearhours  occur  during  the  previous  one. 

The  first  order  statistics  (minimum,  maximum,  mean,  and  standard 

deviation)  were  calculated  for  the  half-hourly  and  the  40HRLP  records 

for  each  variable.  These  are  also  presented  in  the  following  tables. 

In  Section  10  are  plots  of  the  measured  (half-hourly)  temperature  and 

pressure  data.  The  plots  in  Section  11  illustrate  the  detided  pressure 

records,  P.  ,  ,  before  subtracting  the  drift,  with  the  drift  curve 
□etide 

superimposed.  The  P  .  ,  records,  plotted  in  Section  12,  were 

obtained  by  subtracting  the  mean  pressure,  the  drift  curve  (if  nonzero), 

and  the  tide  record  from  the  half-hourly  measured  pressures.  The 

statistics  for  P  ,  ,  are  Included  in  the  following  tables, 

residual 
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TABLE  9.1 

PIES84B2 


lES  Serial  Number:  055 
Pressure  Sensor  Serial  Number:  8181 


Position:  35'’47.81  N 
73*26.99  W 


Depth:  3570  m 


DATE  GMT 

LAUNCH:  Sep  24,  1983  1023 

RECOVERY:  Jan  11,  1984  2004 


CRUISE 

EN106 

0C144 


KEASQRED  PSESStJRE  RECORDS 

(Fig.  10.11) 


DATE 

1st  DATA  POINT:  Sep  24,  1983 

LAST  DATA  POINT:  Nov  18,  1983 


GMT  YEARHOUR 

112952  -2364.5025 

052952  -1050.5025 


Number  of  points:  2629 
Sampling  Interval:  0.50  hrs 


Minimum  =  3623.79  dbar  Mean  =  3624.43  dbar 

Maximum  =  3625.17  dbar  Standard  deviatica  =  0.325  dbar 


RESIDUAL  PRESSURE  RECORDS 


P 


residual 


P 

meas 


MEAN  -  DRIFT  -  TIDE 


DRIFT  =  P,  [1  -  exp(P,  t)]  +  p, 
where  t  =  Time  of  sample  in  hours,  starting  with 
t  =  13.0  for  the  first  data  point. 

Pt  =  -0.255695  dbar 
P,  =  -0.0110545  h  ^ 

P,  =  0.239596  dbar 


TIDE  calculated  from  the  following  constituents: 


M2 

N2 

S2 

K2 

K1 

01 

PI 

01 

H  (dbar): 

.42427 

.10616 

.08304 

.01971 

.09128 

.06666 

.02991 

.01460 

G*: 

353.50 

335.77 

20.90 

21.65 

183.08 

186.63 

182.51 

194.59 

_ DATE  g«I  YEARHOUR 

1st  DATA  POINT:  Sep  24,  1983  232952  -2352.5025 

LAST  DATA  POINT:  Nov  18,  1983  052952  -1050.5025 


Number  of  points:  2605 
Sampling  Interval:  0.50  hrs 


Minimum  =  -0.1059  dbar 
Maximum  =  0.1140  dbar 


Mean  =  0.0000  dbar 
Standard  deviation  =  0.0370  dbar 


TABLE  9.2 
PIES85BCM2 


lES  Serial  Number:  055 
Pressure  Sensor  Serial  Number: 


! 

8181 


Position:  35'’48.09  N 
73'’25.88  W 


Depth:  3560  m 


DATE  GMT 

LAUNCH:  Jan  16,  1984  2344 

RECOVERY:  Jan  17,  1985  0104 


CRUISE 

0C144 

EN124 


MEASURED  PRESSURE  RECORDS 

(Fig.  10.12) 


I 

f 


Minimum 

Maximum 


DATE  GMT 

1st  DATA  POINT:  Jan  17,  1984  005927 

LAST  DATA  POINT:  Jan  16,  1985  235927 


Number  of  points:  17567 
Sampling  Interval:  0.50  hrs 


YEARHOUR 

-8399.009 

383.9915 


=  3645.84  dbar  Mean  =  3646.57  dbar 

=  3647.71  dbar  Standard  deviation  =  0.338  dbar 


RESIDUAL  PRESSURE  RECORDS 


P 


residual 


P 

meas 


MEAN 


DRIFT 


TIDE 


I 


DRIFT  =  P,  [1  -  exp(P,  t)]  +  P, 
where  t  =  Time  of  sample  in  hours,  starting  with 
t  =  13.0  for  the  first  data  point 
P,  =  -0.378324  dbar 
P,  =  -0.0022225  h  ^ 

P,  =  0.359908  dbar 


TIDE  calculated  from  the  following  constituents: 


M2 

N2 

S2 

K2 

K1 

01 

PI 

01 

•w 

.43233 

.10587 

.08715 

.02063 

.09064 

.06984 

.02990 

.01485 

z 

352.84 

333.99 

19.68 

20.286 

181.05 

186.12 

181.76 

184.73 

_ DATE 

1st  DATA  POINT:  Jan  17,  1984 

LAST  DATA  POINT:  Jan  16,  1985 


GMT  YEARHOUR 
125927  -8387.009 
235927  383.991 


Number  of  points:  17543  j 

Sampling  Interval:  0.50  hrs 


Minimum  =  -0.1680  dbar 
Maximum  =  0.1190  dbar 


Mean 

Standard  deviation 


0.0000  dbar 
0.0369  dbar 


! 


Minimum 

Maximum 


Minimum 

Maximum 


40HSLP  PKESSOSE  BECOSDS 

(Fig.  12.5) 


DATE 

1st  DATA  POINT;  Jan  19,  1984 
LAST  DATA  POINT:  Jan  15.  1985 


GMT  YEARHOUR 

000000  -8352.0000 

180000  354.0000 


Number  of  points:  1452 
Sampling  Interval:  6.00  hrs 


=  -0.150  dbar  Mean  =  0.0000  dbar 

=  0.087  dbar  Standard  deviation  =  0.0334  dbar 


TEMPEEATOEE  EECOEDS 

(Figs.  10.2a  and  10.2b) 


DATE 

1st  DATA  POINT:  Jan  17  1984 

LAST  DATA  POINT:  Jan  16,  1985 


GMT  YEARHOUR 
125927  -8387.009 
235927  383.991 


Number  of  points:  17543 
Sampling  Interval:  0.50  hrs 


=  2.166  “C 

=  2.435  °C 


Mean 

Standard  deviation 


2.232  “C 
0.080  "C 


40HELP  TEMPEEATUEE  EECOEDS 

(Fig.  12.2) 


_ DATE 

1st  DATA  POINT:  Jan  x9,  1984 

LAST  DATA  POINT:  Jan  15,  1985 


GMT  YEARHOUR 

000000  -8352.0000 

180000  354.0000 


Number  of  points:  1452 
Sampling  Interval:  6.00  hrs 


Minimum  =  2.168  “C 
Maximum  =  2.432  ®C 


Mean  =  2.234  ®C 
Standard  deviation  =  0.054  °C 


TABLE  9.3 

PIES85BCK3 


Minimum 

Maximum 


lES  Serial  Number:  034 
Pressure  Sensor  Serial  Number:  18426 

Position:  SSMl-OO  N  Depth:  3930  m 

73‘'08.02  W 

DATE  GMT  CRUISE 

LAUNCH:  Jan.  15,  1984  0352  0C144 

RECOVERY:  Washed  up  on  beach  in  Bermuda 


MEASUKED  pbessdse  kecoeds 

(Figs.  10.13a  and  10.13b) 


DATE 

1st  DATA  POINT:  Jan.  15,  1984 

LAST  DATA  POINT:  Jan.  3,  1985 


GMT  YEARHOUR 

04  51  25  -  8443.143 

02  51  25  52.356 


Number  of  points:  16991 
Sampling  Interval:  0.50  hrs 


(Pressure  record  Is  Jumpy.) 


TEMFESATtJKE  KECOKDS 

(Figs.  10.3a  and  10.3b) 


DATE 

1st  DATA  POINT:  Jan.  15,  1984 

LAST  DATA  POINT:  Jan.  3,  1985 


GMT  YEARHOUR 

045125  -  8443.143 

035125  52.356 


Number  of  points:  16991 
Sampling  Interval:  0.50  hrs 


=  2.440  *C  Mean  =  2.467  “C 

=  10.719  ‘C  Standard  deviation  =  0.107  “C 
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TABLE  9.4 
P1ES85C0 


lES  Serial  Number:  053 
Pressure  Sensor  Serial  Number: 


17911 


Position:  36*25.25  N 

73*19.75  W 

DATE 

LAUNCH:  Jan  18.  1985 

RECOVERY:  May  12.  1985 


Depth:  3310  m 


CRUISE 

EN124 

EN130 


MEASIJSED  PSESSUSE  SECOSDS 

(Fig.  10.14) 


1st  DATA  POINT: 
LAST  DATA  POINT: 


DATE 

Jan  18.  1985 
May  12.  1985 


GMT 

210004 

143004 


429.0011 

3158.5011 


Minimum 

Maximum 


3342.76  dbar 
3344.35  dbar 


Number  of  points:  5460 
Sampling  Interval:  0.50  hrs 

Mean 

Standard  deviation 


RESIDUAL  PRESSURE  RECORDS 


3343.51  dbar 
0.335  dbar 


P  ,  .  ,  =  P  -  mean  -  TIDE 

residual  meas 

TIDE  calculated  from  the  following  constituents: 


H  (dbar): 
G*: 


M2  N2  S2  K2 

.43132  .10572  .09210  .02238 

351.03  334.45  18.33  19.56 


1st  DATA  POINT: 
LAST  DATA  POINT: 


Jan  19.  1985 
May  12.  1985 


.  K1  01  PI  01 

3  .09088  .07117  .03012  .01475 

3  182.00  185.80  182.65  183.86 

GMT  YEARHOUR 

090004  441.0011 

143004  3158.0011 


Minimum  = 
Maximum  = 


Number  of  points:  5436 
Sampling  Interval:  0.50  hrs 

-0.1731  dbar  Mean 

0.1222  dbar  Standard  deviation 


0.0000  dbar 
0.0417  dbar 


mikr 


so 
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.j 
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40HBLP  PKESSOSE  BECOBDS 

(Fig.  12.6) 


1st  DATA  POINT: 
LAST  DATA  POINT: 


_ DATE 

Jan  20.  1985 
May  11,  1985 


GMT 

180000 

060000 


474.0000 

3126.0000 


Minimum  = 
Maximum  = 


-0.0802  dbar 
0.0831  dbar 


Number  of  points:  443 
Sampling  Interval:  6.00  hrs 

■  Mean 

■  Standard  deviation 


TEKPEKATOSE  BECOBDS 

(Fig.  10.4) 


0.0000  dbar 
0.0369  dbar 


1st  DATA  POINT: 
LAST  DATA  POINT: 


Jan  18.  1985 
May  12,  1985 


GMT 

210004 

143004 


429.0011 

3158.5011 


Minimum  =  2.300  *C 

Maximum  =  10.436  ®C 


Number  of  points:  5460 
Sampling  Interval:  0.50  hrs 

Mean 

Standard  deviation 


40HBLP  TEXPEBATUBE  BECOBDS 

(Fig.  12.3) 


=  2.373  ®C 
=  0.140  ®C 


1st  DATA  POINT: 
LAST  DATA  POINT: 


Jan  20.  1985 
May  11.  1985 


GMT 

180000 

060000 


474.0000 

3126.0000 


Minimum  =  2 . 312  *C 
Maximum  =  2.461  *C 


Number  of  points:  443 
Sampling  Interval:  6.00  hrs 

Mean 

Standard  deviation 


2.369  “C 
0.032  ®C 


TABLE  9.5 

PIES84C1 


lES  Serial  Number:  056 
Pressure  Sensor  Serial  Number:  17848 


Position:  36“17.20  N 

73*11.40  W 


Depth:  3450  m 


DATE  GMT 

LAUNCH:  Nov  01.  1983  1903 

RECOVERY:  Jan  11.  1984  1459 


CRUISE 

EN107 

0C144 


MEASURED  PRESSURE  RECORDS 

(Fig.  10.15) 


_ DATE 

1st  DATA  POINT:  Nov  01.  1983 

LAST  DATA  POINT:  Jan  11.  1984 


QMT  YEARHOUR 
200406  -1443.9316 
143406  254.5683 


Number  of  points:  3398 
Sampling  Interval:  0.50  hrs 

Minimum  =  3513.51  dbar  Mean  =  3514.20  dbar 

Maximum  =  3515.04  dbar  Standard  deviation  =  0.331  dbar 


RESIDUAL  PRESSURE  RECORDS 


P  ,  =  p  ^  -  mean  -  TIDE 

residual  measured 

TIDE  calculated  from  the  following  constituents: 


M2 

N2 

S2 

K2 

K1 

01 

PI 

01 

H  (dbar): 

.42659 

.09910 

.08669 

.02037 

.09116 

.06876 

.03045 

.01264 

G* : 

353.61 

335.45 

21.82 

23.81 

181.16 

188.76 

182.42 

185.16 

DATE 

gMI.  . 

YEARHOUR 

1st  DATA 

POINT:  Nov 

02.  1983 

080406 

1431.9316 

LAST  DATA 

POINT:  Jan 

11.  1984 

143406 

254.5683 

Number  of  points: 

3374 

Sampling 

Interval: 

0.50  hrs 

Minimum 

=  -0.1433 

dbar 

Mean 

=  0.0000 

dbar 

Maximum 

=  0.1374 

dbar 

Standard  deviation 

=  0.0405 

dbar 

40HSLP  PSESSDBE  BECOBDS 

(Fig.  12.4) 


Minimum 

Maximum 


Minimum 

Maximum 


DATE 

GMT 

YEARHOUR 

1st  DATA  POINT:  Nov 

3.  1983 

180000 

-1398.0000 

LAST  DATA  POINT:  Jan 

10.  1984 

060000 

222.0000 

Number  ' 

Of  points: 

271 

Sampling 

Interval: 

6.00  hrs 

=  -0.0828  dbar 
=  0.0818  dbar 


Mean 

Standard  deviation 


0.0000  dbar 
0.0341  dbar 


XEBPEBATOHE  SECOSDS 

(Fig.  10.5) 


DATE 

1st  DATA  POINT:  Nov  2.  1983 

LAST  DATA  POINT:  Jan  11.  1984 


GMT  YEARHOUR 
080406  -1431.9316 
143406  254.5683 


Number  of  points:  3374 
sampling  Interval:  0.50  hrs 

=  2.219  *0 
=  2.349  “C 


Mean  =  2.257  *C 
Standard  deviation  =  0.038  ®C 


40HSLP  TEXPERATUKE  RECORDS 

(Fig.  12.1) 


_ EAIE _ 

1st  DATA  POINT:  Nov  3.  1983 

LAST  DATA  POINT:  Jan  19.  1984 


GMT  YEARHOUR 

180000  -1398.0000 

060000  222.0000 


Number  of  points;  271 
Sampling  Interval:  6.00  hrs 


Minimum  =  2.221  "C 
Maximum  =  2,347  “C 


Mean  =  2.257  “C 
Standard  deviation  =  0.029  ®C 


TABLE  9.6 

PIBS85CCM1 


lES  Serial  Number:  056 
Pressure  Sensor  Serial  Number:  17848 

Position:  36*15.23  N  Depth:  3475  m 

73*09.89  W 


LAUNCH: 

RECOVERY: 


DATE  GMT  CRUISE 
Jan.  17.  1984  0505  OC144 
Jan.  14.  1985  0120  EN124 


(Pressure  sensor  did  not  work..) 


TABLE  9.7 

PIES85C1 


lES  Serial  Number:  035 
Pressure  Sensor  Serial  Number:  17849 


Position:  36*15.26  N 
73*09.70  W 


Depth:  3475  m 


DATE  GMT 

LAUNCH:  Jan  14.  1985  0217 

RECOVERY:  May  12,  1985  1243 


CRUISE 

EN124 

EN130 


MEASXrSED  PSESSUSE  SECOSDS 

(Fig.  10.16) 


_ date 

1st  DATA  POINT:  Jan  14,  1985 

LAST  DATA  POINT:  May  12.  1985 


GMT  YEARHOUR 

032453  315.4147 

122453  3156.4147 


Number  of  points:  5683 
Sampling  Interval:  0.50  hrs 


Minimum  =  3529.66  dbar  Mean  *  3530.39  dbar 

Maximum  =  3531.21  dbar  Standard  deviation  =  0.333  dbar 


RESIDUAL  PRESSURE  RECORDS 


P  ,  =  P  -  MEAN  -  TIDE 

residual  meas 


TIDE  calculated  from  the  following  constituents: 


M2 

N2 

S2 

K2 

K1 

01 

PI 

01 

.43174 

.10568 

.09249 

.02246 

.09047 

.06960 

.03003 

.01404 

351.20 

334.39 

18.03 

19.09 

182.40 

185.18 

182.91 

183.49 

DATE 

1st  DATA  POINT:  Jan  14.  1985 

LAST  DATA  POINT:  May  12,  1985 


GMT  YEARHOUR 
152453  327.4147 
122453  3156.4147 


Number  of  points:  5659 
Sampling  Interval:  0.50  hrs 


Minimum  =  -0.0891  dbar 
Maximum  =  0.1065  dbar 


Mean  =  0.0000  dbar 
Standard  deviation  =  0.0310  dbar 


40HSLP  PBESSOSE  SECOSDS 

(Fig.  12.6) 


Minimum 

Maximum 


Minimum 

Maximum 


Minimum 

Maximum 


_ PATE 

1st  DATA  POINT:  Jan  16,  1985 

LAST  DATA  POINT:  May  11.  1985 


GMT  YEARHOUR 

000000  360.0000 

060000  3126.0000 


Number  of  points:  462 
Sampling  Interval:  6.00  hrs 


=  -0.0632  dbar 
=  0.0763  dbar 


Mean  =  0.0000  dbar 
Standard  deviation  =  0.0271  dbar 


TEMPERATURE  RECORDS 

(Fig.  10.6) 


_ DATE 

1st  DATA  POINT:  Jan  14,  1985 

LAST  DATA  POINT:  May  12,  1985 


YEARHOUR 

032453  315.4147 

122453  3156.4147 


Number  of  points:  5683 
Sampling  Interval:  0.50  hrs 


=  2.188  "C  Mean  =  2.280  “C 

=  8.291  ®C  Standard  deviation  =  0.111  “c 


40HRLP  TEMPERATURE  RECORDS 

(Fig.  12.3) 


_ PATE 

1st  DATA  POINT:  Jan  16.  1985 

LAST  DATA  POINT:  May  11,  1985 


GMT  YEARHOUR 

000000  360.0000 

060000  3126.0000 


Number  of  points:  462 
Sampling  Interval:  6.00  hrs 


=  2.190  *0  Mean  =  2.277  ®C 

=  2.387  ®c  Standard  deviation  =  0.042  ®C 


TABLE  9.8 

PIES84CC1I2 


lES  Serial  Number:  057 
Pressure  Sensor  Serial  Number:  17849 


Position:  36''05.02  N 

72'’59.94  W 


Depth:  3660  m 


DATE 

LAUNCH;  Nov  1,  1983 

RECOVERY:  June  7,  1984 


GMT  CRUISE 

2158  EN107 

1514  EN118 


MEASURED  PRESSURE  RECORDS 

(Fig.  10.17) 


DATE 

1st  DATA  POINT:  Nov  1.  1983 

LAST  DATA  POINT:  Jun  7,  1984 


GMT  YEARHOUR 

230740  -1440.8723 

150740  3807.1277 


Number  of  points:  10497 
Sampling  Interval:  0.50  hrs 

Minimum  =  3732.73  dbar  Mean  =  3733.56  dbar 

Maximum  =  3734.59  dbar  Standard  deviation  =  0.354  dbar 


RESIDUAL  PRESSURE  RECORDS 


P 


residual 


P  -  MEAN  -  DRIFT  -  TIDE 
meas 


DRIFT  =  P,  [1  -  exp(P,  t)]  +  P, 
where  t  =  Time  of  sample  in  hours,  starting  with 
t  =  13.0  for  the  first  data  point 
P,  =  -0.513066  dbar 
P,  =  -0.0002825  h"'^ 

P,  =  0.246240  dbar 


TIDE  calculated  from  the  following  constituents: 


M2 

N2 

S2 

K2 

K1 

01 

PI  . 

01 

H  (dbar): 

.43285 

.10601 

.08994 

.02138 

.09200 

.06898 

.03032 

.01438 

G® : 

352.23 

332.50 

19.29 

19.71 

180.70 

185.78 

181.45 

183.89 

DATE  GMT  YEARHOUR 

1st  DATA  POINT:  Nov  2,  1983  110740  -1428.8723 

LAST  DATA  POINT:  Jun  7,  1984  150740  3807.1277 


Number  of  points:  10473 
Sampling  Interval;  0.50  hrs 


Minimum  =  -0.190  dbar 
Maximum  =  0.118  dbar 


Mean  =  0.0000  dbar 
Standard  deviation  =  0.0449  dbar 


40HSLP  PBESSDSE  SECOEDS 

(Fig.  12.4  and  12.5) 


1st  DATA  POINT: 
LAST  DATA  POINT: 


3.  1983 
6,  1984 


GMT 

180000 

060000 


-1398.0000 

3774.0000 


Number  of  points:  863 

Sampling  Interval:  6.00  hrs 


Minimum  = 
Maximum  = 


-0.1658  dbar 
0.0916  dbar 


Mean  =  0.0000  dbar 
Standard  deviation  =  0.0416  dbar 


TEMPEBATDSE  EECOEDS 

(Fig.  10.7) 


1st  DATA  POINT: 
LAST  DATA  POINT: 


Nov  2,  1983 
June  7,  1984 


GMT 

110740 

150740 


-1428.8723 

3807.1277 


Number  of  points:  10473 
Sampling  Interval:  0.50  hrs 


Minimum  =  2.244  *C 
Maximum  =  2.482  *C 


Mean  =  2.317  *0 
Standard  deviation  =  0.076  *C 


40HSU>  TEHPEKATDSE  SECOSDS 

(Fig.  12.1  and  12.2) 


1st  DATA  POINT: 
LAST  DATA  POINT: 


Nov  3,  1983 
Jun  6.  1984 


GMT 

180000 

060000 


-1398.0000 

3774.0000 


Minimum  =  2.246  ‘C 
Maximum  =  2.475  *0 


Number  of  points:  863 

Sampling  Interval:  6.00  hrs 

Mean 

Standard  deviation 


2.319  'C 
0.058  'C 


TABLE  9.9 

PIES85CC1I2 


lES  Serial  Number:  054 
Pressure  Sensor  Serial  Number:  8180 


Position:  36*05.07  N 

72*59.86  W 


Depth:  3660  m 


DATE 

LAUNCH:  June  7.  1984 

RECOVERY:  May  12,  1985 


GMT  CRUISE 

1705  EN118 

1004  EN130 


MEASURED  PRESSURE  RECORDS 

(Fig.  10.18) 


DATE 

1st  DATA  POINT:  June  7,  1984 

LAST  DATA  POINT:  May  12.  1985 


GMT  YEARHOUR 

181007  -4973.8314 

094007  3153.6686 


Number  of  points:  16256 
Sampling  Interval:  0.50  hrs 


Minimum  =  3732.42  dbar  Mean  =  3733.34  dbar 

Maximum  =  3734.33  dbar  Standard  deviation  =  0.344  dbar 


RESIDUAL  PRESSURE  RECORDS 


P 


residual 


P 

meas 


MEAN  -  DRIFT  -  TIDE 


DRIFT  =  P,  [1  -  exp(P,  t)]  +  P,  +  P.  t 
where  t  =  Time  of  sample  in  hours,  starting  with 
t  =  13.0  for  the  first  data  point 
Pj  =  0.308608  dbar 

P,  =  -0.0036410  h~^ 

P,  =  -0.407702  dbar 
P,  =  0.000027  dbar/h 


TIDE  calculated  from  the  following  constituents: 


M2 

N2 

S2 

K2 

K1 

01 

PI 

01 

H  (dbar): 

.43341 

.10592 

.08873 

.02120 

.08852 

.06976 

.02935 

.01447 

G*: 

352.60 

335.00 

19.74 

20.77 

181.13 

186.05 

181.82 

184.82 

_ DATE 

GMT 

YEARHOUR 

1st  DATA  POINT:  June  8,  1984  061007  -4961.8314 

LAST  DATA  POINT:  May  12.  1985  094007  3153.6686 


Number  of  points:  16232 
Sampling  Interval:  0.50  hrs 


i 


V 


Minimum  =  -0.1120  dbar 
Maximum  =  0.1215  dbar 


Mean  =  0.0000  dbar 
Standard  deviation  =  0.0350  dbar 
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40HKLP  PEESSOSE  RECORDS 

(Fig.  12.5  and  12.6) 


DATE 


1st  DATA  POINT: 
LAST  DATA  POINT: 


June  9,  1984 
May  11.  1985 


Number  of  points: 
Sampling  Interval: 


GMT 

120000 

000000 

1343 

6.00  hrs 


YEARHOUR 

-4932.0000 

3120.0000 


Minimum  = 
Maximum  = 


-0.082  dbar 
0.099  dbar 


Mean  =  0.0000  dbar 
Standard  deviation  =  0.0315  dbar 


TEMPERATURE  RECORDS 

(Fig.  10.8) 


DATE 


1st  DATA  POINT: 
LAST  DATA  POINT: 


June  7,  1984 
May  12,  1985 


Number  of  points; 
Sampling  Interval: 


GMT 

181007 

094007 


16256 
0.50  hrs 


YEARHOUR 

-4973.8314 

3153.6686 


Minimum  =  2.204 
Maximum  =  4.619 


Mean 

Standard  deviation 


40HRLP  TEMPERATURE  RECORDS 

(Fig.  12.2  and  12.3) 


2.256 

0.071 


DATE 


1st  DATA  POINT: 
LAST  DATA  POINT: 


June  9,  1984 
May  11,  1985 


GMT 

120000 

000000 


YEARHOUR 

-4932.0000 

3120.0000 


Minimum  =  2.205  ®C 
Maximum  =  2.334  ®C 


Number  of  points:  1343 
Sampling  Interval:  6.00  hrs 

Mean 

Standard  deviation 


2.257  ®C 
0.036  ®C 


'.V. 

Sf  J 


TABLE  9.10 

PIES84CCX3 


lES  Serial  Number:  036 
Pressure  Sensor  Serial  Number:  17911 


Position: 

35*48.42  N 
72*42.55  W 

Depth:  3900  m 

DATE 

GMT 

CRUISE 

LAUNCH: 

RECOVERY: 

Jan  15.  1984 
June  7,  1984 

0822 

2109 

0C144 

EN118 

MEASBEED  PSESSDEE  KECOEDS 

(Fig.  10.19) 


DATE 

1st  DATA  POINT:  Jan  15.  1984 

LAST  DATA  POINT:  Jun  7,  1984 


GMT  YEAEHOUR 

092953  345.4981 

205953  3812.9981 


Number  of  points:  6396 

Sampling  Interval:  0.50  hrs 


Minimum  =  3990.15  dbar  Mean  =  3990.94  dbar 

Maximum  =  3991.77  dbar  Standard  deviation  =  0.337  dbar 


EESIDUAL  PEESSUEE  EECOEDS 


P 


residual 


P 


measured 


MEAN 


TIDE 


TIDE  calculated  from  the  following  constituents: 


M2 

N2_ 

S2 

K2 

K1 

01 

PI 

01 

H  ( dbar ) : 

,43048 

.10519 

.09131 

.02181 

,09100 

.06813 

.02987 

.01475 

G*: 

352.05 

332.17 

19.27 

19.72 

181.50 

185.05 

181.98 

184.06 

_ date _  GMT  YEARHOUR 

1st  DATA  POINT:  Jan  15,  1984  212953  357.4981 

LAST  DATA  POINT:  Jun  7,  1984  205953  3812.9981 


Number  of  points:  6912 

Sampling  Interval:  0.50  hrs 


Minimum  =  -0.1641  dbar  Mean  =  0.0000  dbar 

Maximum  =  0.1060  dbar  Standard  deviation  =  0.0529  dbar 


40HSLP  PEESSTJEE  EECOEOS 

(Fig.  12.4  and  12.5) 


Minimum 

Maximum 


Minimum 

Maximum 


_ DATE 

1st  DATA  POINT:  Jan  17,  1984 

LAST  DATA  POINT:  Jun  6.  1984 


GMT  YEARHOUR 

060000  390.0000 

120000  3780.0000 


Number  of  points:  566 

Sampling  Interval:  6.00  hrs 


=  -0.1363  dbar 
=  0.0881  dbar 


Mean  =  0.0000  dbar 
Standard  deviation  =  0.0511  dbar 


TEKPEEATDEE  EECOEDS 

(Fig.  10.9) 


_ DATE 

1st  DATA  POINT:  Jan  15,  1984 

LAST  DATA  POINT:  Jun  7,  1984 


GMT  YEARHOUR 
212953  357.4981 
205953  3812.9981 


Number  of  points:  6912 

Sampling  Interval:  0.50  hrs 


=  2,365  *C 
=  2.494  "C 


Mean  =  2.396  ®C 
Standard  deviation  =  0.039  'C 


40HELP  TEHPEEATUSE  EECOEDS 

(Fig.  12.1  and  12.2) 


_ DATE 

1st  DATA  POINT:  Jan  17,  1984 

LAST  DATA  POINT:  Jun  6,  1984 


GMT  YEARHOUR 

060000  390.0000 

120000  3780.0000 


Number  of  points:  566 

sampling  Interval:  6.00  hrs 


Minimum  =  2.367  ®C 
Maximum  =  2.486  “C 


Mean  =  2.397  “c 
Standard  deviation  =  0.019  °C 


TABLE  9.11 
PIES85CCM3 


lES  Serial  Number:  058 
Pressure  Sensor  Serial  Number:  19327 

Position:  35*48.23  N  Depth:  3890  m 

72*42.57  W 


DATE 

GMT 

CRUISE 

LAUNCH:  June  7,  1984 

RECOVERY:  May  12,  1985 

2239 

0635 

EN118 

EN130 

MEASURED  PRESSURE  RECORDS 

(Fig.  10.20) 

DATE 

GMT 

YEARHOUR 

1st  DATA  POINT:  June  7,  1984 

LAST  DATA  POINT:  May  12,  1985 

235920 

062920 

-4968.0111 

3150.4889 

Number  of  points: 
Sampling  Interval: 

16238 

0.50  hrs 

Minimum  =  3988.67  dbar  Mean  =  3989.38  dbar 

Maximum  =  3990.26  dbar  Standard  deviation  =  0.334  dbar 


RESIDUAL  PRESSURE  RECORDS 


P 


residual 


P  -  MEAN  -  DRIFT  -  TIDE 
meas 


DRIFT  =  P,  [1  -  exp(P,  t)]  +  P, 
where  t  =  Time  of  sample  in  hours,  starting  with 
t  =  13.0  for  the  first  data  point 
P,  =-0.196539  dbar 
P,  =-0.0012165  h  ^ 

Pj  =  0.177363  dbar 


TIDE  calculated  from  the  following  constituents: 


M2 

N2 

S2 

K2 

K1 

01 

PI 

01 

H  (dbar): 

.43211 

.10560 

.08885 

.02127 

.08816 

.06896 

.02922 

.01424 

G“: 

352.83 

335.33 

20.12 

21.21 

181.19 

186.55 

181.94 

185.23 

:lSi 


m 


I  » 


DATE 


1st  DATA  POINT: 
LAST  DATA  POINT: 


June  8, 
May  12, 


1984 

1985 


Number  of  points: 
Sampling  Interval: 


115920 

062920 


16214 
0.50  hrs 


YEARHOUR 

-4956.0111 

3150.4889 


Minimum  =  -0.121  dbar 
Maximum  =  0.148  dbar 


Mean 

Standard  deviation 


0.0000  dbar 
0.0369  dbar 


'Si 

§ 

‘o' 

S 


40HKLP  PSESSnSE  BECOBDS 

(Fig.  12.5  and  12.6) 


_ DATE 

1st  DATA  POINT:  June  9.  1984 

LAST  DATA  POINT:  May  11.  1985 


_2M3L_  YEARHOUR 
180000  -4926.0000 
000000  3120.0000 


Number  of  points:  1342 
Sampling  Interval:  6.00  hrs 


Minimum  -  -0.100  dbar  Mean  =  0.0000  dbar 

Maximum  =  0.120  dbar  Standard  deviation  =  0.0334  dbar 


TEHPEBATDSE  RECORDS 

(Fig.  10.10) 


1st  DATA  POINT: 
LAST  DATA  POINT: 


Minimum  =  2.375  *C 
Maximum  =  5.983  *C 


_ DATE 

June  7,  1984 
May  12.  1985 

Number  of  points: 
Sampling  Interval: 


-gtOL  YEARHOUR 
235920  -4968.0111 
062920  3150.4889 

16238 
0.50  hrs 


Mean  =  2.414  ®C 
Standard  deviation  =  0.077  *C 


40HRLP  TEXPERATURE  RECORDS 

(Fig.  12.2  and  12.3) 


_ PATE 

1st  DATA  POINT:  June  9.  1984 

LAST  DATA  POINT:  May  11.  1985 


gtq.  YEARHOUR 

180000  -4926.0000 

000000  3120.0000 


Number  of  points:  1342 
Sampling  Interval:  6.00  hrs 

Minimum  =  2.377  "c 
Maximum  =  2.512  *C 


Mean  =  2.415  “C 
Standard  deviation  =  0.024  °C 
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10.  HALF-HOURLY  DATA  FOR  EACH  INSTRUMENT 


j 

Plots  of  the  measured  bottom  pressure  and  temperature  are 
presented.  The  time  scale  is  the  same  for  all  plots,  with  each 

increment  corresponding  to  5  days.  The  axis  begins  on  0000  GMT  of  the  | 

first  date  labelled. 

The  vertical  scale  for  each  variable  is  consistent  between 

instruments.  Each  increment  corresponds  to  0.5  dbar  for  the  bottom  i 

o 

pressure  data  and  to  0.02  C  for  the  temperatures. 

The  sampling  interval  is  0.5  hours  for  all  instruments.  The  length 
and  start  and  end  times  of  the  data  records  are  tabulated  in  the  i 

previous  section. 

10-1 _ Temperature  Data 

The  raw  temperature  records  are  shown  in  Figures  10.1  to  10.10. 

All  the  temperatures  are  in  the  range  2.2  to  2.5°C,  with  a  tendency  for 
the  offshore  (deeper)  records  to  be  warmer.  We  have  observed  this 
pattern  in  previous  investigations  in  this  region,  where  the  colder  core 
of  the  Western  Boundary  Undercurrent  (WBUC)  is  somewhat  onshore  of 
(shallower  than)  the  array.  The  temperature  records  are  characterized 
by  "square-pulse-liRe"  warming  events;  i.e.,  they  have  many  plateaus  of 


relatively  constant  temperature  (±0.02°C)  separated  by  very  rapid 
changes  (-0.05  to  0.15°C  within  a  few  hours)  to  a  different  plateau  that 
may  last  typically  for  1  to  10  days.  We  tentatively  associate  this 
behavior  with  lateral  shifts  of  the  offshore  temperature  front  of  the 
WBUC,  although  it  may  also  be  caused  by  advectlve  pulses  of  warmer  or 


colder  water. 


10.2  Pressure  Data 

The  raw  pressure  records  are  shown  in  Figures  10.11  to  10.20.  The 
semidiurnal  tidal  oscillations  of  roughly  1  dbar  range  are  prominent  in 
these  records,  and  the  29-day  repeat-cycle  of  spring  and  neap  tides  is 
readily  observed.  Other  features  of  these  pressure  records,  which  are 
only  slightly  noticeable  on  this  coarse  plotting  scale,  are  discussed  in 


Section  12,  after  the  tides  and  the  drift  have  been  removed  from  them. 


Figure  10.2b  Kaw  temperature  record  of  PIES85BCM2,  continued 


bBCM3 


Raw  temperature  record  of  PIES85BCM3,  continued 
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pressure  record  of  PIES85BCM2 


Figure  10.16  Raw  pressure  record  of  PIES85C1 
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11.  DETIDED  PRESSURE  RECORDS  AND  DRIFT  CURVES 


£ 


% 


The  detided  pressure  (P,  .  )  records,  which  have  had  the  mean  and 

detide 

the  tides  removed,  are  presented  in  Figures  11.1  to  11.10.  For  the 
records  which  had  a  measurable  long-term  drift,  the  regressed  drift 
curves  are  superimposed  on  the  pressure  curves. 

The  vertical  scale  of  0.1  dbar  per  increment  is  consistent  between 
all  records,  except  for  PIES85BCM3  (Figure  11.3),  which  had  large 
unexplained  jumps.  For  that  instrument,  the  vertical  scale  has  been 
doubled  to  0.2  dbar  per  increment. 

The  time  scale  is  the  same  for  all  plots  with  each  increment 
corresponding  to  25  days.  The  axis  begins  on  0000  GMT  of  the  first  date 
labelled. 

The  sampling  Interval  is  0.5  hours  for  all  instruments.  The  length 
and  start  and  end  times  of  the  data  are  tabulated  in  Section  9. 
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Figure  11.5  Detlded  pressure  record  of  PIES84C1 
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Figure  11.7  Detlded  pretsur*  record  of  PIES84CCM2 


PRESSURE  (DBAR) 


12.  FULLY  PROCESSED  RECORDS 


The  low-pass  filtered  temperature  and  bottom  pressure  records  are 
presented  for  each  instrument.  The  40HRLP  temperature  records  are  shown 
first.  These  are  followed  by  the  40HRLP,  detided,  dedrifted  pressure 
records . 

The  time  scale  is  the  same  for  all  plots,  with  each  increment 
corresponding  to  10  days.  The  axis  begins  on  0000  GMT  of  the  first  date 
labelled.  Since  the  measurements  at  several  of  the  sites  were  obtained 
by  two  or  more  instruments,  the  data  for  January  to  May  1984  is  repeated 
in  the  figures  to  show  coinciding  measurements.  Each  record  is  labelled 
with  the  Instrument  names  above  the  appropriate  portions. 

The  vertical  scale  for  each  variable  is  consistent  between 

instruments.  Each  increment  corresponds  to  0.05  dbar  for  the  bottom 

o 

pressure  measurements  and  0.04  C  for  the  temperatures. 

The  sampling  interval  is  6  hours  for  all  the  records.  The  length 
and  start  and  end  times  of  the  data  records  are  tabulated  in  Section  9. 
12.1  Temperature.:.  ^OHRLP  data 

The  40HRLP  temperature  records  are  shown  in  Figures  12.1  to  12.3. 
the  square-pulse  nature  of  the  temperature  records  has  been  rounded 
somewhat  by  the  low-pass  filtering.  Visually,  the  temperature  records 
from  different  sites  have  overall  only  a  small  amount  of  coherence; 
many  pulses  are  quite  different  from  one  site  to  another. 


coalescence  into  the  Gulf  Stream  (Bane  and  Watts,  1985).  Pressure  and 
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the  first  six  days  of  deployment.  ^ 

We  estimate  the  residual  error  in  the  final  pressure  records,  after  the  "dedrifting" 
calculations,  to  be  tv-pically  0.02  dbar  r.m.s.  {or  0.06  dbar  r.m.s.)  if  the  first  6 
days  of  the  record  are  excluded  (or  included,  respectively) .  This  low  drift-error 
opens  many  possibilities  for  studies  that  require  knowledge  of  the  low-frequency 
dynamic  pressure  signal  in  the  deep  ocean. 

Part  I  on  Methods  contains  a  short  review  of  bottom  pressure  measurement  in  the  deep 
ocean,  a  description  of  the  sensors  that  we  used,  a  discussion  of  their  performance 
and  drift  relative  to  t>'pe  of  construction  and  prior  pressurization  history  ("precon¬ 
ditioning"),  and  estimates  of  the  accuracy  of  the  dedrifted  pressure  records. 

In  Part  II  of  this  report,  the  full  data  processing  is  described,  including  calibration 
parameters,  corrections  for  the  influence  of  temperature  variations  on  the  pressure 
sensor,  and  parameterization  to  remove  sensor  drift  errors  by  least-squares  regression 
onto  an  exponentially  decaying  time-dependence.  Time  series  are  plotted  which 
illustrate  several  steps  in  the  processing:  the  edited  half-hourly  pressure  records, 
the  detided  pressures  with  drift-model  curves  superimposed,  and  the  low-pass  filtered, 
"dedrifted"  pressure  records  (i.e.,  after  subtracting  the  estimated  drift  curve). 


